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In a recent Letter [l[ Efetov et al. propose an ex- 
act mapping of an interacting fermion system onto a 
new model that is supposed to allow sign-problem free 
Monte Carlo simulations. In this Comment, we show 
that their formalism is equivalent to the standard ap- 
proach of Blackenbecker, Scalapino and Sugar (BSS) Q 
for fermionic systems and has the same sign statistics and 
minus sign problem. 

Our first observation is that the partition function for 
a given configuration of the auxiliary fields <f> is the same 
in the standard formulation Zf [Eq. (8) in Ref. [l[] and 
in their new bosonized scheme Z b [Eq. (9)]: 


z f [</>] = z b [</>] . 


(i) 


This observation is trivial in the limit of the time step 
At — > 0, where both schemes reproduce the same par- 
tition function. Since Zf can be negative also in this 
limit 0, Zf, is also not sign-positive. Both Zf and Z b 
are positive if 4> is a smooth path Q , but restricting the 
configuration space to smooth paths amounts to a semi- 
classical approximation. 

We next show that Zf [</>] = Z b [<fi] also holds for finite 
At and piecewise constant paths where the field 4> r (T) 
is constant on the interval [(/ + 1)A, /A]. Efetov et al.'s 
Eq. (9) is equivalent to Z b [<j)] = Tr (e-^°[//(/3, 0)) 
■£-Ui(n,T ) = - #7,1 (n)Z7i (ti,7-o) and 
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We are left with 
th 


the task of evaluating Ui (/A, (/ — 1)A) on the I time 
interval where the fields are constant and take the 
values (j) r i. For this constant in time field configura- 
tion, Ui(lA, (I - 1)A) = e *AiTo e -Aft !e -(;-l)Affo with 

Hi = H - J2 r <fir,im z , r . Hence, 


Z b [(/)} = det 


N 


i + n< 


-Ah, 


1=0 


(2) 


We also demonstrate that the proposed formalism is 
equivalent to that of BSS Starting from the expres- 
sion for Z found between Eq. (9) and Eq. (10) and using 
a matrix notation the partition function reads 


In 


Zb[<t>] 
Z n 


du^Tr<$>G a 


(3) 


where $ is a matrix with elements $(r/,r'i') = 
4>r,i5 rr '8w , the trace is over the space time indices and 


G a = G° + uaG°<S>G a = (l- uaG°<P) 1 G°, 


(4) 


where G° is the non interacting Green function. From 
Eqs. © and @ we get 


Z J Q 4" 1- M <tG°$ 


(5) 


The integral in terms of u and a can be carried out ana- 
lytically yielding the BSS partition function 


Z b [0] = Z exp (-Tr [in (l - (G°$) 2 )] ) 

= -- = n 1 , 

det(l - (G°$) 2 ) 11 detG CT 


(6) 


A sign problem is present in Eq. ([3]), but is hidden in 
the u integral. If G°$ has eigenvalues on the real axis 
with absolute value greater than 1, then the u integral 
runs over poles. This will give phases of iir in the expo- 
nent which can lead to negative values of the exponential. 

Finally, Efetov's et al. suggest an alternative way 
of evaluating the same weight Z b from a bosonic field 
A rr / (t); however, the equation of motion that A satisfies, 
[their Eq. (13)], is singular. Furthermore, the equation 
of motion cannot uniquely determine A — not even with 
their additional constraint J2 r -^ rr ( r ) = 0- The latter 
can be shown, for example, using a two site cluster and 
a single time slice. 

In conclusion, we have shown that for piecewise contin- 
uous paths the partition function obtained from Efetov 
et al.'s method [1| is equivalent to the standard BSS for- 
mulation and has, in particular, the same sign problem. 
However, their method provides a new perspective on 
the minus sign problem, as it can be viewed as originat- 
ing from poles of the coupling constant (u) integral, or 
as branch cuts of a logarithm when the integral over u is 
performed analytically. 
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